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Abstract 

Spin-dependent tunneling through an indirect bandgap barrier Uke the GaAs/AlAs/GaAs het- 
erostructure along [001] direction is studied by the tight-binding method. The tunneUng is char- 
acterized by the proportionality of the Dresselhaus Hamiltonians at T and X points in the barrier 
and by Fano resonances (i.e. pairs of resonances and anti-resonances or zeroes in transmission). 
The present results suggest that large spin polarization can be obtained for energy windows that 
exceed significantly the spin splitting. The width of these energy windows are mainly determined 
by the energy difference between the resonance and its associated zero, which in turn, increases 
with the decrease of barrier transmissibility at direct tunneling. 

We formulate two conditions that are necessary for the existence of energy windows with large 
polarization: First, the resonances must be well separated such that their corresponding zeroes are 
not pushed away from the real axis by mutual interaction. Second, the relative energy order of the 
resonances in the two spin channels must be the same as the order of their corresponding zeroes. 

The degree to which the first condition is fulfilled is determined by the barrier width and the 
longitudinal effective mass at X point. In contrast, the second condition can be satisfied by choosing 
an appropriate combination of spin splitting strength at X point and transmissibility through the 
direct barrier. 

PACS numbers: 72.25-b,72.10.Bg,73.21.Fg,73.21.Ac,73.23.-b 



I. INTRODUCTION 



The spin rather than the charge of carriers has attracted a lot of interest leading to a 
new field of electronics dubbed spintronics.— In this context, spin-polarized transport in 
non-magnetic semiconductor structures and spin-dependent properties originating from the 
spin-orbit interaction are a promising road to spin based devices.— 

Despite the progress that has been made,'^ spin injection from ferromagnetic leads proved 
to be very challenging.— Consequently, spin-dependent transport in nanostructures comprised 
of non-magnetic semiconductors has been the focus of extensive work in the past years.— 
Recent theoretical research has suggested that the current resulting from electron tunneling 
through zinc-blende semiconductor single^i^ or double barrier structures^ can be highly spin 
polarized. The origin of the spin-dependent tunneling in these structures stems from the 
fact that the barrier material lacks center of inversion. 

In the effective mass approximation, the electron Hamiltonian of a zinc-blende structure 
has an additional spin-dependent coupling called the Dresselhaus term^^ 



where cr, are the Pauli matrices, and k^, ky, and kz are the components of electron wave 
vector. For a barrier along [001] direction the Dresselhaus Hamiltonian is reduced to 



Perel, Tarasenko, and coworkers^i^ showed that the F point Hamiltonian in Eq. Q induces 
an effective mass correction leading to a spin-polarized transmission. It is important to 
emphasize that the spin dependent part of the effective mass Hamiltonian at X point isiSiii 



and therefore proportional to the Hamiltonian in Eq. (j2)). The Hamiltonians in (j2|) and (jH} 
are diagonahzed by 



Hd = 1 [(t.K [kl - kl) + Gyky [kl - kl) + a^kz [kl - kl)] , 
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FIG. 1: Tunneling through GaAs/ AlAs/ GaAs barrier. T band edges are shown by solid Une and 
X band edges are represented by dashed Une. Ei and E2 are X-valley quasi-bound states in the 
AlAs barrier. Direct tunneling (r-F-r) is shown by solid curly arrow and F -X-T tunneling is 
depicted by dotted curly arrow. 

with the polar angle of the wave vector kt in xy plane. Therefore the spin states are not 
mixed by the interaction between X states and F states in the barrier. 

Spin-dependent transport can be studied using numerous treatments such as the k ■ p 
approach, full-band tight-binding calculations, and ab initio methods. For various reasons 
the theoretical study of spin tunneling through an indirect barrier like GaAs/ AlAs/ GaAs 
has not been fully addressed before, k-p cannot fully address the problem because the AlAs 
barrier accommodates at least one quasi-bound state into X-valley. Thus, beside the F-F-F 
tunneling, which occurs through the higher F-valley, one must also consider the tunneling 
(F-X-F) through the lower X-barrier (see Fig. P). The k ■ p method is a perturbative 
method that can be "tuned" for the necessities of spin-dependent processes (for instance, see 
Ref. Q, in which spin- dependent evanescent states in the band gap are studied). In contrast, 
empirical tight-binding methods provide a treatment of the full Brillouin zone, but they lack 
the complete description of the Dresselhaus term when the spin-orbit is introduced.-^ This 
is due to the fact that the orthogonality assumption in tight binding models is incompatible 
with the formulation of the spin-orbit interaction.^^ In principle, the above shortcomings 
should be overcome by utilizing ab-initio density functional theory methods. However, 
these methods suffer on the side of bandgap reproducibility.-^ 

Spin dependent tunneling has been recently analyzed with a 1-band envelope-function 
model.-^ In their study,-^ the authors neglected the (F-F-F) tunneling and the presence of 
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X-valley quasi-bound states in the AlAs barrier. However, spin tunneling through the indi- 
rect barrier of the GaAs/AlAs/GaAs heterostructure shows another pecuhar property. The 
confined X states in the AlAs barrier interact^"^ with the continuum F states in GaAs forming 
Fano resonancesiS (i.e.:pairs of resonances/anti-resonances). In this paper we demonstrate 
that one can use the proximity in the resonance and anti-resonance states in conjunction 
with the spin splitting produced by the spin-dependent Hamiltonian to obtain a large degree 
of spin polarization within the range between the resonance and anti-resonance energy. For 
this purpose we devise a spin-dependent tight-binding model that provides a realistic view 
of the spin-dependent tunneling through an indirect barrier. We convert the spin- dependent 
effective mass Hamiltonians for a single band to their tight-binding versions following the 
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recipes of Ref. Q. The coupling between F and X valleys is made according to Ref. 

The paper is organized as follows. In next section a simple model is analyzed in order 
to gain insight into the physics of spin- dependent tunneling. The third section contains a 
realistic tight-binding model and the numerical results. Conclusions are drawn in the fourth 
section. 



II. TIGHT-BINDING MODEL OF SPIN-DEPENDENT TUNNELING 

Consider a simple tight-binding (TB) model of the spin- dependent tunneling through an 
indirect barrier. The main assumptions for this TB model are: spin states are degenerate 
in left/right lead (bulk-like states) and the Dresselhaus Hamiltonians are proportional for F 
and X states in the barrier, so we can assume that the spin states in the leads are eigenvalues 
of the Dresselhaus Hamiltonian. 

a. Spin tunneling through an indirect barrier. The Hamiltonian of the system is 

-1 

H = E [£C^,aCn,a + {t C+_^^„Cn,a + h.C.)] + 

n=—co,a='\,i 
oo 

E [^<aCn,a + {t C+„Cn+l,a + h.C.)]+ , . 

£icJ^Co,t + £^2cJ|Co,i + (VicJ^c_i,-f + V'acJ^c.i,^ + h.c.) + 
(ViC({^^ci,| + V2C^^ci,| + h.c.) + {ticti^^ci^-\ + t2cti^^ci^i + h.c.) 
The first two terms on the right hand side of Eq. (0) are the Hamiltonians of the contacts 
(leads), where e and t are the on-site energy and transfer integral, respectively (spin de- 
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generate ). c^^jcna is the creation/annihilation of an electron with spin a on site n. The 
remaining part is the Hamiltonian of the barrier and its couphng to the leads. The active 
region is modeled by three sites: site —1 that is like the left hand side contact, site that 
is the actual barrier, and site 1 that is like the right hand side contact. Thus the effective 
left /right hand side contact ends/starts at site —2/2. The matrix form of the Hamiltonian 
for the sites -1, 0, and 1 (in fact E — H, where E is energy) with appropriate boundary 
conditions for an open system is^i 
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where S/, /j (E) are the self-energies of the semi- infinite parts, i.e., 

where Hl h are the Hamiltonians of the semi-infinite left and right hand sides and 6 is an 
infinitesimal positive number. The retarded Green function 

of the left/right hand side semi- infinite contact in Eq. ((Tj) is actually the diagonal part 
Gf (E)_2 _2 /Gr (E)^ 2 representing the sites -2/2. The expressions of these Green function 
elements can be found from their equation of motion and the use of the finite difference 
equation method.— If we consider the parameterization e — E = 2t cos (ka), with k a 
complex parameter and a a lattice constant parameter, one obtains the following equation 
for self-energies, 

^L,R = -t e""'. (9) 
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Therefore, the Green function for the sites -1, 0, and 1 with the boundary conditions for 
an open system is 



{E) 



-, -1 



t e-'^" 










-h 








-t e-''"' 





-V2 







-V* 





E-Ei 





-V{ 








-v; 





E-e2 





-v; 


^1 










-t e-'^"" 








''2 





-V2 








(10) 



We notice that the Hamiltonian is not hermitian due to open boundary conditions. To 
calculate the transmission probability from site 1 to we use the formula^! 



T {E, h) = Tl {E, h) iGf^ (E, kt)Y Tr (E, k 



t) , 



(11) 



with 



Tl,r {E, k) = I [Si,^ (E, h) - SI ,j (E, h)] . (12) 

Since the Dresselhaus Hamiltonians are proportional at F and X points in the barrier^iSiii 
we can solve separately for each spin. The Green function for 'spin up' is 

-t e-'''^ -Vi -ti 
Gf{E)= E-ei -V{ . (13) 

-t* -Vi -t e-''"' 

A similar equation is obtained for the 'spin down'. The poles of are the solutions of the 
determinant equation 
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TABLE I: Matrix elements in eV of the nearest neighbor model outlined in Eq. ©. 



Indirect barrier RTD-like structure 



t 1.0 1.0 

Vi 0.05 0.005 

V2 0.052 0.0052 

h 0.05 0.0 

t2 0.052 0.0 

£1 0.17 0.17 

£2 0.175 0.175 



while the zeroes of the transmission are the zeroes of the Green function [G^) ^ ^ relating 
the sites -1 and 1, 



-V^ E-e, 
-ti -V^ 



A 



0. 



(15) 



The equation for poles reads 



(16) 



and for zeroes 



{E-ei)h = -\Vi\\ (17) 

Since t is basically the conduction band bandwidth and ti and Vi are tunneling rates, then 
Vi,ti « t, such that the pole is given by 



The equation for the zero is 



W P 

E^6i- 2^^e''"'. (18) 



\V |2 

E = e,-^-^. (19) 
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Since we have Vi,ti << t, the energy separation between the pole and the zero is about 
Vi/ti. The resonances and zeroes occur at shghtly different energies for the two spin chan- 
nels, resulting in a large spin polarization due to the combination of a sharp increase in 
transmission at resonance followed by an abrupt decrease to zero at anti-resonance. Hence, 
the energy range of large spin polarization will depend on Vi and ti but not on the mag- 
nitude of spin splitting. One can notice that the energy separation between resonance and 
anti- resonance can be increased by decreasing ti, i.e., increasing the width of the barrier. 
In Fig. 121 we illustrate the above arguments with the parameters given in Table □ We also 
calculate the spin polarization with the equation 

Tt - Tl 

p = ^r^- (20) 

For comparison we also plot the case of the resonance tunneling diode (RTD) configuration. 
The RTD configuration is made by setting ti and )f:2 to 0. To separate the spin resonances, 
we also have chosen the values of Vi and V2 ten times smaller than their values in the 
indirect-barrier configuration. Fig. El shows that the indirect barrier configuration has a 
clear advantage over the RTD configuration. 

b. Spin tunneling through a direct barrier. Following the same procedure one can also 
calculate spin transmission through a direct barrier. The Green function of the open system 
is 
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The Green function for the 'spin up' is 
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Using Eqs. (fTT|) and (fT^ we calculate the spin- dependent transmission probabilities for a 
direct barrier. The transmission probabilities are approximated by 
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4|ti| 
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FIG. 2: Spin dependent transmission (upper panel) and spin polarization (lower panel) through 
an indirect barrier versus a resonant tunneling diode configuration. One can see the wide energy 
window of large polarization created by indirect tunneling. 

and 

n-^^sin^ka, (24) 

i.e., the transmission for 'spin up' is different from transmission for 'spin down' leading to 
spin polarization. 



III. RESULTS OF A REALISTIC TIGHT-BINDING MODEL AND DISCUS- 
SIONS 

In order to make quantitative assessments of the spin tunneling through an indirect 
barrier, we devise a spin-dependent tight-binding model for a system that is partitioned into 
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layers. The layers from — oo to and from + 1 to oo are the contacts, while the layers 
from 1 to are the active layers. First, we consider the spin independent Hamiltonian. The 
coupling between T and X states is treated similarly to Ref. |20|, with the Hamiltonian 



H 



(25) 



HxT Hx 

Hy and Hx are the Hamiltonians at F and X point, respectively. H^x and Hxv are the 
couplings between F and X at the interface layers. For simplicity, we do not distinguish 
between Xi and X3, such that Hy and Hx are single-band effective mass Hamiltonians that 
are converted to TB Hamiltonians according to the parameterization given in Ref. 13 This 
TB parameterization has been successfully used in quantum transport for non-equilibrium 
conditions and incoherent scattering processes.^ The parameterizationi^jS^ is made for the 
effective mass Hamiltonian 

where m\ is the effective mass in the left contact, the effective mass is considered z- 
dependent, and the spatial dependence of the transverse energy has been incorporated in 
the transverse momentum (/c^) dependent potential: 



The corresponding tight-binding parameters for the non-diagonal part are: 



^ (m, + m,) A2 ^^^^ 



and the diagonal part is 



A ih) = 7^J— + —)+V^ ih) . (29) 

In Eqs. ()28p and ()29|). rrii is the effective mass at site i on the mesh of spacing A, Vi (kt) is the 
potential at site i, which also includes the band offsets, m~ = ^ and = !Ih+^h+i_ 

The spin dependent Hamiltonian is expressed in the basis spanned by spinors such that 
the Hamiltonian is diagonal in this basis. At F point the spin dependent part of the effective 
mass Hamiltonians is introduced through the corrections to the effective masses defined in 
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TABLE II: Dresselhaus coefficients at F and X points for GaAs and AlAs, calculated with the 
GW method^^'^^ with spin-orbit coupling included. The units are the atomic units. 
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FIG. 3: (a) Spin dependent transmission probability of 1-band (direct) tunneling through 
GaAs/ (AlAs) 8 /GaAs barrier, (b) Comparison between full (direct and indirect) and 1-band (direct) 
transmission probability of a GaAs/ (AlAs)8 /GaAs barrier, (c) Comparison of spin polarization be- 
tween full (direct and indirect) and 1-band (direct) electron transmission of a GaAs/(AlAs)8/GaAs 
barrier, (d) Spin polarization obtained from 1-band calculation over a broader energy range. 

Eq. Q and the effective potential defined in Eq. (jTTj) . The spin-dependent part at X is 
expressed through different band offsets for the two spin projections as one can see from 
Eq. ©. 

Our calculations are performed with the effective masses and band edges taken from 
Ref . 13 The Dresselhaus parameters given in Table |n] are calculated with a quasi-particle 
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FIG. 4: Spin dependent transmission probability and spin polarization of a GaAs/(AlAs)Ar/GaAs 
heterostructure. A^=4, 6, and 8. 

self-consistent GW (G'=Green function, VF=screened Coulomb interaction) method as in 



Ref. 



25 



26l with spin-orbit coupling included perturbatively. The splitting at T in GaAs 
found in the present work with the GW method is three times smaller than the value used 
by Perel et al. in Ref. The GW spin splitting at F for AlAs is about 2.5 times smaller 
than in GaAs. While for GaAs one can also use the experimental estimation for the splitting 
at the F poinlj^l there are no such estimations for AlAs. Therefore to be consistent, we used 
the GW values for both. The value of spin splitting at X point is almost the same as the one 
obtained within the local density approximation (LDA) with our FP-LMTO (full potential- 
linear muffin-tin orbitals) code.^^ The X-F band offset for GaAs/AlAs heterostructure is 
chosen to be 160 meV.-^ Throughout the paper we have chosen a value = ~ ^ 0.05 for 
the transverse wave vector (a is the lattice constant of GaAs). 

In Fig. El we compare the one-band model (direct tunneling) given by the effective mass 



12 



Hamiltonian and the two-band model (direct and indirect tunneling) given by Eq. ()25|) . Be- 
low the indirect barrier, the main contribution to spin tunneling and polarization is provided 
by direct tunneling. However, the tunneling and the polarization through X states become 
dominant for energies above the indirect barrier. The result shows that the confinement of 
the X states in the barrier increases the energy threshold at which the tunneling through X 
states becomes dominant. In the calculations^^ in which the confinement of the X states in 
the barrier are neglected, the tunneling is predominantly indirect for energies slightly below 
the top of indirect barrier. Therefore, our calculations suggest that multi-band calculations 
are needed to fully describe the electron transport in these heterostructures. 

In Fig. 0] we show the transmission probability and spin polarization for 
GaAs/(AlAs)Ar/GaAs heterostructures with = 4, 6, and 8. Energy windows with large 
polarization can be seen between the resonance and its corresponding zero. The width of the 
window increases with the barrier width as it has been demonstrated in the previous section. 
Only the first resonance has a corresponding zero on the real axis, for the other resonances, 
the zeroes are pushed off the real axis.-^ Therefore, if the resonances are close enough, no 
well defined window with large spin polarization can be found. The possibility to obtain 
well separated resonances with zeroes on the real axis is controlled by the combination of the 
longitudinal effective mass in the barrier at X point and the barrier width. A lighter longi- 
tudinal effective mass and/or a narrower barrier push farther apart the resonance energies 
in the barrier. 

In Fig. we plot the transmission probability and spin polarization for 
GaAs/(AlAs)Ar/GaAs with = 5 and 7. The parity of the number of AlAs monolayers 
has been taken into account.— The case A^= 5 shows two wide windows with large and 
opposite spin polarizations. Again, the polarization windows are mainly determined by the 
resonance and anti-resonance positions and not by the magnitude of spin splitting. However, 
N = 7 shows no such energy windows because the zeroes have moved away from the real 
axis. Moreover, at larger values of A^ no energy windows with large polarization are found 
for both even and odd values of A^. 

The pattern can be re-established by increasing the strength of the Dresselhaus coeffi- 
cient at X point in the barrier. A barrier made of Alo.8Gao.2As can achieve this goal. In 
Alo.8Gao.2As, the Dresselhaus coefficients are mixtures of those of AlAs and GaAs. GaAs 
has larger coefficients, in particular, j3 is ten times larger than the (3 coefficient of AlAs, 
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FIG. 5: Spin dependent transmission probability and spin polarization of a GaAs/(AlAs)Ar/GaAs 
heterostructure. N = 5 and 7. 

making the f3 coefficient of the compound stronger than that of AlAs. In Fig. IHl we make 
a comparison between Alo.8Gao.2As and AlAs barriers with a thickness of N=12 monolay- 
ers. Virtual crystal approximation was employed to calculate the physical parameters of 
Alo.8Gao.2As. This is a reasonable assumption, since AlAs and GaAs have similar structural 
and electronic properties. Fig. IHl illustrates clearly that Alo.8Gao.2As barrier shows a window 
of polarization, while AlAs barrier does not. 

Visual analysis of Fig. El and Fig. [7| suggests that in order to obtain polarization windows, 
the resonances and zeroes in the spin channels must be properly ordered. If a resonance 
is occurring ffist in the spin channel 1, the zero in the transmission coefficient of the spin 
channel 1 must precede the zero in the transmission coefficient of spin channel 2. This 
condition is fulffiled for wider barriers provided that the Dresselhaus coefficient at X point 
is sufficiently large. 

The practical aspect of focusing the electrons to energies within the large-polarization 
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FIG. 6: Comparison of spin dependent transmission probability and spin polarization between 
GaAs/(AlAs)i2/GaAs and GaAs/(Alo.8Gao.2As)i2/GaAs heterostructures. 

window can be achieved by placing an RTD structure in front of the tunneling barrier. In 
this way, one not only can control the incoming energy of the electrons but also ensures that 
most electrons have non-vanishing transverse momenta'^- and, consequently, their energies 
are spin-split. 

It is important to analyze at this point the influence that the neglect of the Dresselhaus 
term in the leads can have on the conclusions drawn in the present work. Since the tight- 
binding basis is localized, different basis can be used to treat the spin Hamiltonian in barrier 
and in the leads. In the leads, the basis is that which makes diagonal the Hamiltonian in 
Eq. Q, while in the barrier the basis is given in Eq. (j^. The non-diagonal spin dependent 
part is then transferred to interface terms between leads and barrier. The corrections to the 
Green function and transmission coefficients are quadratic with respect the strength of this 
non-diagonal term. This analysis indicates that the conclusions of the article are not likely 
to change if the Dresselhaus term in the leads is properly taken into account. A detailed 
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FIG. 7: Closer look at spin dependent transmission probability around resonance energy for 
GaAs/(AlAs)i2/GaAs and GaAs/(Alo.8Gao.2As)i2/GaAs heterostructures. It explains the origin 
of energy window with large spin polarization. 

quantitative discussion will be presented elsewhere since special care must be exercised due 
to the presence of the fc^-linear term in Eq. 

IV. CONCLUSIONS 

We investigated the spin dependent transport across an indirect semiconductor barrier of 
a zinc-blende structure like GaAs/AlAs/GaAs heterostructure along [001] axis by means of a 
combination of several tight-binding models. Spin tunneling through such an indirect barrier 
exhibits two major characteristics: the proportionality of the Dresselhaus Hamiltonians at 
r and X points and the Fano resonances. A generic tight-binding Hamiltonian has shown 
that large spin polarization occurs in the energy window determined by the separation 
between the resonance and its associated anti-resonance and not by the magnitude of the 
spin splitting of resonances. Moreover, the energy separation between the resonance and its 
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corresponding anti-resonance increases as the barrier width increases. 

Reahstic calculations have been performed with a two-band tight-binding model. The 
effective mass Hamiltonians at F and X have been convertedi^ to the tight-binding Hamil- 
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tonians. The T-X coupling was implemented following the scheme presented in Ref. 
Accordingly, the Dresselhaus Hamiltonians at F and X in the barrier have been included in 
the effective masses and band offsets. The calculations show that, in order to obtain energy 
windows with large polarization, two conditions need to be satisfied. The first condition 
consists of having well separated resonances such that their corresponding anti-resonances 
do not interact with each other. The second condition is that the relative energy order of 
the resonances in the two spin channels must be the same as the order of their correspond- 
ing zeroes. The first condition is achieved by an appropriate combination of barrier width 
and longitudinal effective mass at X point, while the second condition is accomplished by 
a combination of spin splitting strength at X point and transmissibility through the direct 
barrier. 

Electrons can be focused in the required energy windows by passing them through a 
resonant tunneling diode structure situated in front of the indirect barrier. Using such an 
experimental setup, one could obtain large spin polarization following the procedure of Perel 
et al^ and Glazov et al^. 



Acknowledgments 

This work has been supported in part by NSERC grants no. 311791-05 and 315160-05. 
The authors wish to acknowledge generous support in the form of computer resources from 
the Reseau Quebecois de Calcul de Haute Performance. 



* Electronic address: Ititus . sandu@umontreal . cal 

^ S. Wolf, D. Awschalom, R. Buhrman, J. Daughton, S. von Molnar, M. Roukes, A. Chtchelka- 

nova, and D. Treger, Science 294, 1488 (2001). 
2 I. Zutic, J. Fabian, and S. D. Sarma, Rev. Mod. Phys. 76, 323 (2004). 

^ D. Awschalom, D. Loss, and N. Samarth, eds.. Semiconductor Spintronics and Quantum Com- 
putation (Springer, Berlin, 2002). 

17 



4 E. 1. Rashba, Phys. Rev. B 62, R16267 (2000). 

^ A. Hanbicki, B. Jonker, G. Itskos, G. Kioseoglou, and A. Petrou, Appl. Phys. Lett. 80, 1240 

(2002). 

^ G. Schmidt, D. Ferrand, L. W. Molenkamp, A. T. FiUp, and B. J. van Wees, Phys. Rev. B 62, 
R4790 (2000). 

^ V. I. Perel, S. A. Tarasenko, I. N. Yassievich, S. D. Ganichev, V. V. Belkov, and W. Prettl, 

Phys. Rev. B 67, 201304(R) (2003). 
8 S. A. Tarasenko, V. I. Perel, and I. N. Yassievich, Phys. Rev. Lett. 93, 056601 (2004). 
^ M. M. Glazov, P. S. Alcksccv, M. A. Odnoblyudov, V. M. Chistyakov, S. A. Tarasenko, and 

I. N. Yassievich, Phys. Rev. B 71, 155313 (2005). 
° G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

^ E. L. Ivchenko and G. E. Pikus, Superlattices and Other Heterostructures. Symmetry and Optical 

Phenomena (Springer, BerHn, 1997), 2nd ed. 
^ N. Rougemaille, H.-J. Drouhin, S. Richard, G. Fishman, and A. K. Schmid, Phys. Rev. Lett. 

95, 186406 (2005). 
3 T. B. Boykin, Phys. Rev. B 57, 1620 (1998). 

T. Sandu, Phys. Rev. B 72, 125105 (2005). 
^ R. Jones and O. Gunnarsson, Rev. Mod. Phys. 61, 689 (1989). 
6 S. Mishra, S. Thulasi, and S. Satpathy, Phys. Rev. B 72, 195347 (2005). 
^ R. Beresford, L. F. Luo, W. I. Wang, and E. E. Mendez, Appl. Phys. Lett. 55, 1555 (1989). 
^ U. Fano, Phys. Rev. 124, 1866 (1961). 

^ W. R. Frensley, in Heterostructures and Quantum Devices, edited by W. Frensley and N. Ein- 

spruch (Academic Press, San Diego, CA, 1994), p. 273. 
'° Y. Fu, M. Willander, E. L. Ivchenko, and A. A. Kiselev, Phys. Rev. B 47, 13498 (1993). 

S. Datta, Electronic Transport in Mesoscopic Systems (Cambridge University Press, New York, 

1995). 

'2 A. B. Chen, Y. M. Lai-Hsu, and W. Chen, Phys. Rev. B 39, 923 (1989). 
'3 R. Lake, G. Khmeck, R. C. Bowen, and D. Jovanovic, J. Appl. Phys. 81, 7845 (1997). 
'"^ G. Klimeck, R. C. Bowen, T. B. Boykin, and T. A. Cwik, Superlatt. and Microstruct. 27, 519 
(2000). 

S. V. Faleev, M. van Schilfgaarde, and T. Kotani, Phys. Rev. Lett. 93, 126406 (2004). 

18 



A. N. Chantis, M. van Schilfgaarde, and T. Kotani (2005), cond-mat/0508274. 
G. E. Pikus, V. A. Marushchak, and A. N. Titkov, Sov. Phys. Semicond. 22, 115 (1988). 
M. Methfessel, M. van Schilfgaarde, and R. Casali, in Electronic Structure and Physical Proper- 
ties of Solids: The Uses of the LMTO Method, Lecture Notes in Physics, edited by H. Dreysse 
(Springer- Verlag, Berlin, 2000), p. 114. 

E. E. Mendez, E. Calleja, and W. I. Wang, Appl. Phys. Lett. 53, 977 (1988). 

R. C. Bowen, W. R. Prensley, G. Klimeck, and R. Lake, Phys. Rev. B 52, 2754 (1995). 

T. Sandu, G. Klimeck, and W. P. Kirk, Phys. Rev. B 68, 115320 (2003). 



19 



